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Abstract 

The meson screening mass in the pseudoscalar channel is calculated from the mo- 
mentum dependent meson spectral function, using HTL fermionic propagators. A 
careful subtraction procedure is required to get an UV finite result. It is shown that 
in the whole range of temperatures explored here the HTL screening mass stays 



above the non- interacting result, slowly approaching the value Tn scr — 2y / 7T 2 T 2 + m^, 
where m 2 ^ is the HTL asymptotic thermal quark mass. Our analysis leads to a bet- 
ter understanding of the excitations of QGP at sufficiently large temperatures and 
may be of relevance for interpreting lattice results. 
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1 Introduction 



This work is devoted to the study of mesonic screening masses, which can be 
denned as the inverse screening length characterizing the exponential falloff 
of the mesonic spatial correlator. The importance of their evaluation in order 
to identify the relevant degrees of freedom in hot and dense QCD was first 
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pointed out in [1,2,3], where the screening masses of mesons and nucleons were 
calculated on the lattice. More recent results can be found in [4,5,6,7,8,9]. 

Though obtained in different discretization schemes, some common features 
characterize the results of all the above studies. At high temperature chiral 
simmetry appears to be restored: n and a, p and a±, N + (even-parity nucleon) 
and N_ (odd-parity nucleon) turn out to be degenerate. As the temperature 
increases meson and baryon screening masses approach their ideal gas value 
of 2irT and 3irT, respectively. Finally, while the screening mass of the vector 
mesons approaches the non-interacting result immediately above the phase 
transition (hence, the presence of resonant states in this channel seems to be 
ruled out), the result for the 7r and a mesons stays below the ideal gas value 
up to higher temperatures!"*! 

Thus, 7r and a states could possibly survive in the deconfmed phase as col- 
lective excitations and be identified with the soft modes of the chiral phase 
transition, as first proposed in [10] . In this connection, the presence of meson- 
like excitations above T c has been also investigated within an effective PNJL 
model [11]. 

Which are then the active degrees of freedom in the QGP phase, leading to the 
correct interpretation of the data obtained both in the heavy-ion experiments 
and in lattice simulations? 

For what concerns energy densities not too far from the deconfinement transi- 
tion, the ones presently achievable in the heavy-ion collision experiments, the 
possible presence of a huge set of neutral and colored bound states surviving 
up to temperatures of order ~ 2T C , has been proposed [12,13,14,15] as be- 
ing able to give a unified interpretation of different experimental evidences: 
large cross sections, small mean free paths, perfect fluid behavior, fast ther- 
malization, elliptic flow, etc. Actually neither the theoretical description of 
the QGP thermodynamics in terms of bound states [16,17], nor the interpre- 
tation of the experimental data in terms of fast thermalization and perfect 
fluid behavior [18] is universally accepted. 

While the features of the system just above T c are still under debate, at higher 
temperatures, starting from T ~ 3T C , the relevant degrees of freedom should 
be weakly interacting quarks and gluons. In this regime, resummation schemes 
based on the Hard Thermal Loop (HTL) approximation were developed. Dif- 
ferent thermodynamical observables turned out to be well described [19,20,21] 
(in particular the slow approach of the entropy, pressure and energy density 
to their ideal gas value is well reproduced) in terms of resummed quark and 
gluon propagators. The latter imply a rich structure of many-body phenom- 



1 Indeed in [9] it was pointed out that, in order to get reliable results for pion 
properties on the lattice, one has to employ a fermionic action which preserves 
chiral simmetry. 
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ena: thermal masses, collective excitations (plasminos and longitudinal gluons) 
and Landau damping. 

It is then of interest to investigate what these dressed degrees of freedom, 
which, for large enough temperatures, nicely reproduce the QGP thermody- 
namics, imply for the correlation of hadronic excitations. 

In this paper, extending the work done in [24,25] where the HTL meson spec- 
tral functions (MSFs) at zero and finite momentum were evaluated, we address 
the calculation of the z— axis correlator of mesonic currents. As we will show, 
from the large distance behavior of the correlator one can extract the screen- 
ing mass of the meson. We limit ourselves to the pseudoscalar channel where, 
as in [24,25], the computation "simply" requires the convolution of two HTL 
resummed quark propagators. 

Though the evaluation of hadronic screening masses is a major achievement 
of lattice studies of the degrees of freedom characterizing the hot QCD (the 
large number of lattice sites available along the spatial directions allows one 
to study the large distance behavior of the correlators), analytical approaches 
are not so common in the literature. 

The large distance behavior of hadronic correlators in the high temperature 
limit was first discussed in [26,27]. The full analytical study of the spatial 
mesonic correlations was first addressed (to our knowledge) in the non in- 
teracting case in Ref. [28]. Effects of the interaction on the meson screening 
masses were analyzed, for example, in Refs. [29,30,31,32,33] within a dimen- 
sional reduction framework. A positive correction of order g 2 T to the non- 
interacting result = 2nT was found within this kind of approach. We will 
compare our numerical results to the above mentioned ones. 
A study of the screening masses in the N JL model was presented in [34] , whose 
main result (beyond stressing the huge numerical difficulty of the calculation) 
is to show how in the scalar and pseudoscalar channels the value of the screen- 
ing mass in the deconfined phase stays below the non-interacting result, thus 
reflecting the presence of non-trivial correlations. This seems to agree with the 
lattice results and to support the presence of pion and sigma excitations for 
temperatures slightly larger than T c , representing the soft modes of the chiral 
phase transition. 



Our paper is organized as follows. We shortly review in Sec. 2 the basic defi- 
nitions of the mesonic operators, correlators and spectral functions addressed 
in this work. In Sec. 3 we discuss how to get an UV finite result for the z-axis 
correlator in the non-interacting case and then we present the results obtained 
from the finite momentum HTL meson spectral function. Finally, in Sec. 4, 
we summarize our results and discuss how they compare to the ones obtained 
in independent approaches, both in the continuum and on the lattice. 
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2 Mesonic spatial correlation function 



An interesting quantity to extract informations on the properties of hot QCD 
is the correlator of currents carrying the proper quantum numbers to cre- 
ate/destroy mesons. In practice, in the imaginary time formalism (which allows 
to study equilibrium properties in thermal field theory), one has to compute 
the following thermal expectation value 

GmHt, x) = (J M (-ir, x)Jl(0, 0)) , (1) 

which provides information on how fluctuations of mesonic currents are cor- 
related. In the above r G [0,(3= 1/^1 > while 

J M (-ir, x) = J M (-ir, x) - (J M (-ir, x)) , (2) 

denotes the fluctuation of the current operator 

J M (-ir, x) = q(-ir, x)T M q(-ir, x) , (3) 

being T M = 1, 7 5 , 7 M , 7^7 5 for the scalar, pseudoscalar, vector and pseudovec- 
tor channels, respectively. 

The correlator in Eq. (1) is conveniently expressed through its Fourier com- 
ponents according to: 

G M {-vr,x) = - £ —^e-^e^ x G M (icu n ,p), (4) 

P n=-oo J O") 

where cu n = 2nnT (n = 0, ±1, ±2 . . . ) are the bosonic Matsubara frequencies. 

The correlator in Fourier space can be written in term of its spectral density 
through the following representation 

+ °° a (lj p) 1 
G M {ioJn,p) = - / duj— — - — =>• ff M (w,p) = -ImG M (w + jf?,p), (5) 

J lUl n — LJ 7T 

— oo 

showing the link between the Meson Spectral Function (MSF) a M and the 
retarded correlator. 

On general grounds one expects the large distance correlations to be expo- 
nentially suppressed, the coefficient governing this damping being identified 
with the mass of the (mesonic) excitation, not necessarily corresponding to 
a bound state. For example the asymptotic high temperature value 27rT, al- 
ready mentioned in the Introduction for the meson screening mass, refers to 
the propagation of a non-interacting qq pair, each of the two particles carrying 
the minimal fermionic Matsubara frequency ttT. 
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Before proceeding in the calculations we wish to stress that the introduction 
of a thermal bath establishes a privileged reference frame. Lorentz invariance 
is then lost and the result obtained from the correlation along the r-axis (usu- 
ally referred to as dynamical mass) will not necessarily be equal to the one 
arising from the z-axis correlator (screening mass). The latter is a quantity 
which is most easily extracted from lattice calculations, through the analysis 
of the exponential damping of the following correlator: 



Gm{z) = J dr J dxs_G M (-iT,Xj_,z) , (6) 
o 

which for large distances is expected to display a behavior like: 

G M (z) ~ e~ m ^ z . (7) 

z — >+oo 

We notice that the integrations in Eq. (6) select the components in Fourier 
space corresponding to vanishing Matsubara frequency and transverse mo- 
mentum. Hence one can write: 



+ 00 

Qu(z)= J =^&* z G„(po = 0,p±=0,p,) 

— oo 
+oo 

I dP^ e ip z z f + °° (L; <7m{u,P±=0,Pz 

J 2w J-oo uj 



where in the second line use has been made of the spectral representation 
given in Eq. (5). This will be the starting point for the present investigation 
of z-axis correlations. 



The high-energy behavior of the spectral function, <Tm{uj) ~ u 2 , for uj — > oo 
makes the integration over u in Eq. (8) UV divergent. This problem is already 
present at the level of the non-interacting theory and has to be cured through 
a proper subtraction procedure. In the next section we show how we have 
solved this problem when quark and antiquarks propagating in the QGP are 
described by HTL propagators, o 



3 Mesonic screening masses: HTL results 



Before addressing the HTL calculation of the z-axis meson correlator, we 
briefly recall the essential aspects of the non-interacting result presented in [28]. 

2 The first complete calculation in which this difficulty was overcome in the non- 
interacting case can be found in [28]. 
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The free correlator is written as the sum of two terms: 



s l0) (z) = s&(z) + s!£L(z)- (9) 

The first one, temperature independent, describes the process in which the 
external probe with the quantum numbers of a meson excites a qq pair from the 
vacuum. The second one accounts for the presence of a medium (the thermal 
bath) which, on the one side, introduces a Pauli-blocking factor for the qq 
excitation; on the other side, it allows a second process (not possible in the 
vacuum), namely the absorption of the external probe by a (anti-) quark from 
the thermal bath, which is then promoted to an unoccupied single-particle 
state. 

The expression for the vacuum piece (employing Pauli-Villars regularization 
in the intermediate steps) is found to be (in the case Nf = 2): 

S&(z) = ^m*K 2 (2m z), (10) 

where mo is the current mass of the quark and K 2 is the modified Bessel 
function of order 2. The matter part, in turn, reads: 





-moK^moz) + ttT ]T exp f-2z ^{21 + 1)% 2 T 2 + m§ 

l=-oo ^ 



, (11) 



where K\ is the modified Bessel function of order 1. 

By summing the two contributions an exact cancellation of the vacuum part 
with the first term of the matter one occurs: notably, being mo <C T, this term 
represents the major contribution to <? ma tt at large distances. One then finds: 

gl°\z) = ^ (1 - - z ^j J>xp (-2*^(2/ + l)W* + mg) , (12) 

which for large distances displays the following asymptotic behavior: 



g(o) (z) „ 2 J^Lj^ T 2 + m 2 e -2^ + ml (13) 

2 ^oo nz v 

This leads to the non-interacting result for the meson screening mass, which 
reads: 

mS = + ml (14) 

and coincides with 2nT for vanishing current quark masses. This result arises 
from a dramatic cancellation of two terms, which otherwise would completely 
dominate the correlator at large distances. The possibility of analytically per- 
forming all the calculations in the free case allows to extract in a clean way 
the exact large distance behavior of the correlator. 
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In the interacting case the problem is much more involved, since the meson 
spectral function is obtained only numerically. /,From Eq. (8) it appears that 
the z-axis correlator can be expressed in terms of the finite momentum MSF. 
In [25] we evaluated the latter in the HTL approximation in the pseudoscalar 
channel, which amounts to a convolution of two HTL resummed quark propa- 
gators. In the following we will limit ourselves to this case. Several non-trivial 
many-body processes turn out to contribute to the MSF and we refer the 
reader to ref. [25] for a detailed discussion of the various terms. 
We then write the interacting spectral function as follows: 



HTL 



[U>,P Z )=0T(<*,P,) - (°5?>,P.) - ^M L (^)) 

= o^{w,p z )-a**(w,p z ), (15) 



where we introduced the auxiliary MSF <r^ x . For the latter we made the fol- 
lowing choice. We employ the expression one gets in the free case from the con- 
volution of two massive quark propagators, leading for the z-axis correlator to 
the result given in Eq. (12). For the mass of the quarks we replace the free value 
m by the HTL asymptotic thermal quark mass = 
m q being the thermal gap mass appearing in the fermionic spectrum. Indeed, 
in the HTL approximation, the asymptotic thermal mass governs the large- 
momentum regime of the quark dispersion relation. Our choice guarantees a 
convergent high-energy behavior for the difference o"^ ff (o;,p z ), which makes 
the integration of the latter in Eq. (8) well defined. 

In Figs. 1-3, we show the u dependence of all the contributions to Eq. 15 for 
a few momenta p z , for the case T = 2T C . The numerical calculations were 
performed up to cu = 2500 fm _1 in order to check the correct high-energy limit 
of the difference between the HTL and the auxiliary spectral functions. Such 
a difference is well behaved and, as u grows, goes smoothly to zero, as we can 
see in Figs. 4-6. 

The z-axis correlator accordingly reads: 

g UTL (z) = g anx (z)-g diS (z), (16) 

where Q anx (z) is given by Eq. (12) after making the replacement m — > m^, 
while Q dlS (z) is obtained by numerically evaluating 

+oo +oo 

G dm (z)= J ^e^ z G^( Pz )= J f^^/J^4f(^), (17) 

— oo — oo 

the integration over u being now well defined. 

/From Eq. (17) it clearly appears that the full determination of Q dlS (z) would 
in principle require performing the Fourier Transform (FT) of a function that 
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Fig. 1. The finite momentum pseudo-scalar MSF divided by u together with the 
auxiliary one and their difference, for p z = 0.6 fm -1 and T = 2T C . 
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Fig. 2. The finite momentum pseudoscalar MSF divided by lo together with the 
auxiliary one and their difference, for p z = 4.0 fm -1 and T = 2T C . 
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Fig. 3. The finite momentum pseudo-scalar MSF divided by oj together with the 
auxiliary one and their difference, for p z = 50.0 fm _1 and T = 2T C . 




Fig. 4. The difference between the auxiliary and the HTL pseudoscalar MSF divided 
by to for momenta p z = 20 — 100 fm _1 and T = 4T C . 
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Fig. 6. The same as in Fig 5 but for large to. 
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we know only numerically, after summing a huge set of different contribution 
which are listed in Ref. [25]. In the present work, however, we are only in- 
terested in the determination of the coefficient which governs the exponential 
damping of Q dlS (z). Hence we will limit ourselves to evaluate with high preci- 
sion the large z behaviour of the latter. 

For this purpose it appears then more convenient to employ the following pro- 
cedure: we perform a fit of G dlS (p z ) with a function whose FT in coordinate 
space is known and displays an exponential decay at large distances. For the 
sake of simplicity, for the fitting function we choose 

G dis ( Pz ) = -j:^- 2 , as) 

which leads to a z-axis correlator with the following functional form: 

Q^{z) = -Y j c l e~ m ^ . (19) 
i=i 

The above, indeed, is a simple multi-exponential fit, which is generally used 
to extract the masses of physical states with given quantum numbers from the 
study of lattice correlators. 

We are aware that there is a priori no reason, for the above expression, to 
be the correct functional form of G dlS (z), as it can be guessed, for example, 
by looking at the free result in Eq. (12). This assumption could obviously 
introduce some bias in our results. Hence, while we are confident that at 
large distances the correlator is dominated by a single decreasing exponential, 
controlled by the smallest between the two parameters rrii, we do not attach 
a real physical meaning to the other coefficients of the fitting function. 

In order to clarify our fitting procedure let us start by considering the large 
distance behavior of the non-interacting z-axis correlator. iFrom Eq. (12), for 
the case m = 0, it follows that the latter can be conveniently expressed in 
the following form: 



(z) = Cl (z)e^ T > + c 2 (z)e-^ z + ... 

= Cl (z)e'^ T >(l + ^rle-^ T > + ...) , (20) 

V °i\ z ) J 

where the Ci(z) are functions displaying only a mild dependence on z (com- 
pared to an exponential), which in principle can be obtained by matching the 
above expression with Eq. (12). Then one expects the z-dependence being 
dominated by the first exponential for distances such that 

(4ttT)^»1 <=► z*>-±=. (21) 
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2 


13.415 


0.015 


30-^100 


4 


26.213 


0.043 


40^-200 


10 


65.84 


0.41 


200^300 



Table 1 

The coefficient mi extracted from the fit of G dlS (p z ) for the different temperatures 
explored, together with the range of momenta over which the fit has been performed. 

We believe it to be a reasonable assumption that these considerations remain 
essentially true also in the interacting case. Then, from 

+oo 

g^ {z) = [ ?El jPzZ G ^( pz) (22) 
J 2n 

— oo 

it follows that in order to correctly reproduce the large z behavior of Q (z) one 
needs a careful determination of G(p z ) up to p z ~ 1/z. Contributions from 
p z ^> 1/z are in fact suppressed by the oscillations of the integrand. Hence, 
from Eq. (21), one can assume that the numerical evaluation of the screening 
mass from a fit reconstruction of G(p z ) requires a sufficiently large number of 
data points for p z < AnT . On the other hand one should not extend the fit to 
too large values of p z , otherwise the lighest mass would no longer control the 
^-behavior. 

Of course a more precise determination of the range over which the fit has to be 
performed cannot be obtained through these semi-quantitative considerations. 
In practice, in fitting G dlS (p z ) we adopted the following procedure. For each 
temperature we varied the maximum value of p z in a "reasonable range" . After 
checking that the value of the lightest mass (mi) was sufficiently robust with 
respect to the arbitrary choice of p max (i.e. displaying very modest variations), 
for each temperature we took a weighted average of the different values. 

The results for mi are reported in Table 1, together with the estimate of the 
error resulting from the fitting procedure and the range over which we varied 
pmax a furt^gr check of the robustness of the values of m x thus obtained, 
for each temperature, we extended our fit to larger values of p z , introducing a 
third decreasing exponential. In all the cases we found that the values of mi 
obtained with this three-mass fit are indeed consistent with the ones given in 
Table 1. 

The accuracy of the fit is shown in Figs. 7-10 for different temperatures. In all 
figures the numerical calculation of G dlS (p z ) (open circles) is compared with 
the smooth fitting function. 

In Table 2 we collect our results, at different temperatures, for the screen- 



12 



-0.14 



-0.16 



-0.18 



a. 




-0.22 



-0.24 



-0.26 



-0.28 



Fig. 7. G dm (p z ) [in fur 2 ] for the temperature T/T c = 1. 



-0.35 



-0.45 



C5 




-0.55 



-0.65 



30 40 
Pz (fm" 1 ) 



80 
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13 




120 



Fig. 9. G dlS (p z ) [in fm 2 ] for the temperature T/T c = 4. 
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Fig. 10. G diS (p z ) [in fm" 2 ] for the temperature T/T c = 10. 
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Table 2 

The screening masses governing the exponential decay of the non-interacting cor- 
relator, of the auxiliary one, of the difference between the HTL and the auxiliary 
correlators, respectively. In the last column we show the HTL screening mass. 



T/T c 


ci ± 5c\ (fm 3 ) 


iri2 i 5ni2 (fm 1 ) 


c 2 ± 5c 2 (far 3 ) 


1 


0.364±0.001 


195.5±4.0 


15.3±0.3 


2 


1.754±0.003 


291.3±1.7 


54.4±0.3 
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10.13±0.03 


377.6±1.7 


208.3±0.8 


10 


126.1±1.5 


527.4±6.6 


1301.0±10.3 



Table 3 

The weighted average of the other fit coefficients extracted from G dlS (p z ) for the 
different temperatures explored. 

ing mass: together with the non-interacting result (mf ree = 2ttT) we report 
the auxiliary HTL mass m™ = 2y / 7r 2 T 2 + m 2 ^ and the lightest mass aris- 
ing from our fit of Q dlS (z). Clearly the screening mass (m^ L ), governing the 
exponential decay of the HTL z-axis correlator G nTL (z), is given by: 

m^ = min(rrC',m 1 ). (23) 

The Table shows that at all temperatures m\ turns out to be smaller than 
but larger than mf ree ; hence it has to be assumed as the "true" HTL 
screening mass. Notice that, independently on the result of the fit, rri^ sets 
an upper cutoff to the HTL screening mass. 

For the sake of completeness, we report in Table 3 the average values found 
for the other fit parameters, though we do not attach any physical meaning to 
them, in particular to the heavier mass m,2- Notice that, since the parameter 
m 2 is always very large, by limiting the range of the fit to not too large values 
of p z one can indeed fit equally well the data by a constant (which would 
correspond to a contact term in coordinate space) plus a single decreasing 
exponential. We checked that in this case the value of m 1 results slightly 
larger than the one we quote in Tables 1 and 2 arising from the two- mass fit, 
and approach m^f even faster. 

In Fig. 11 we plot the HTL and auxiliary screening masses in units of the non- 
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Fig. 11. The temperature dependence of the HTL result for the meson screening 
mass, together with the auxiliary one, in units of the non-interacting value 2irT. 

interacting value 2-rrT as a function of the temperature. One can see that for 
large temperatures, the HTL screening mass approaches the auxiliary result 



m scr ~m scr 



in 



oo ' 



(24) 



where is the asymptotic thermal mass of the quark. The high tempera- 
ture behavior we have found deserves a discussion on its physical origin. The 
main effect of the interaction is to dress the quarks with a thermal mass. Yet 
we remind that the HTL propagator respects chiral invariance, though be- 
ing dominated for large momenta by a quasiparticle pole with the dispersion 
relation of a massive particle. 

It is indeed interesting to evaluate the lowest order correction to the non- 
interacting result for the screening mass, which can be obtained by expanding 
Eq. (24) for <C nT (as it is indeed the case, since = gT/V3). One 
gets, for large temperatures: 



HTL ^ 
"^scr — 



2ttT + — g 2 T . 
3n 



(25) 



One can compare this result with the one obtained in [32,33] which can be 
cast in the form: 

m scr ~2vrT+ (J- + A^j g 2 T , (26) 
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A being an additional correction term, to the same order. 

Our result has been obtained through the convolution of two HTL quark prop- 
agators, somehow assuming that the relevant degrees of freedom in the high- 
temperature phase of QCD are dressed quasiparticles. This leads (assuming 
the asymptotic high-temperature behavior numerically found to be correct) 
to a non-perturbative result for the meson screening mass, which implicitly 
contains infinite powers of g 2 . 

In [32,33] the approach is completely different. These authors do not use re- 
summed propagators. An effective (2 + 1) dimensional Lagrangian is built, 
in which the relevant degrees of freedom are the lowest fermionic Matsubara 
modes and soft gluons. Within this approach the result given in Eq. (26) is 
found for the screening mass, where the coefficient of the g 2 T term results 
from the sum of two pieces. The first one arises from a correction to the mass 
of the fermionic modes and coincides with our result, to the extent that the 
HTL result is expanded perturbatively; the second term instead (A) is related 
to vertex and quark self-energy corrections, arising from soft gluons, which are 
ignored in our approach. Indeed we remind the reader that the pseudoscalar 
vertex receives no HTL correction and, in the evaluation of HTL self-energy 
diagrams, one considers the loop integrals being dominated by hard internal 
momenta. 

Both our approach and the one employed in [32,33] agree in predicting a 
small positive correction of order g 2 T to the non-interacting result for the 
meson screening mass, in contrast with the lattice results available so far. The 
above discussed origin of the small discrepancy, between the two approxima- 
tion schemes, which affects the exact numerical coefficient of g 2 T, makes, in 
our opinion, this qualitative trend of the results rather robust and represents 
a challenge for lattice studies. 

Finally, in Fig. 12 we display the behavior of the spatial correlator Q(z) for the 
case T = 2T C obtained by combining togheter the auxiliary correlator Q a,uyi {z) 
(evaluated from Eqs. (16) and (12) ) and the result for Q dlS (z) arising from 
our two-mass fit. 



4 Conclusions 

In this paper we have addressed the numerical evaluation of the screening mass 
of a (pseudoscalar) meson in the QGP phase of QCD, i.e. the quantity which 
governs the large-distance exponential decay of the correlations of mesonic 
current operators. 
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Fig. 12. The final result for the HTL z-axis correlator Q nTL (z) (continuous line) 
obtained, according to Eq. (16), through the difference between the auxiliary cor- 
relator aux (z) and Q dlS (z), which arises from our simple exponential fit. Only the 
large distance behavior should be trusted. The figure refers to the case T = 2T C 

We have presented results obtained in the Hard Thermal Loop approxima- 
tion. While the pseudoscalar vertex does not receive HTL corrections, for the 
fermionic lines the HTL resummed propagators have been employed. In this 
scheme the spectral density of the quark correlator is characterized by two 
stable quasi-particle poles (normal quark and plasmino mode) for time-like 
momenta, and by a continuum contribution in the space-like domain related 
to the damping of the quark due to its interaction with the other particles of 
the thermal bath (Landau damping). In a previous publication [25] we have 
shown in detail that this structure of the fermionic propagator gives rise to an 
extremely rich and complex set of many-body processes contributing to the 
finite-momentum spectral function of a meson in the QCD deconfined phase. 

Here we employed the previous results to derive the spatial meson correlator, 
limiting the discussion to the asymptotic propagation along the z-axis and 
presenting numerical results for the screening mass of a pseudoscalar meson. 
We discussed in the text how we managed to deal with the problem of the UV 
divergences (arising from the u 2 growth of the MSF at high energy) which 
one encounters in the calculation: the latter are overcome through a careful 
subtraction procedure which takes advantage of how this problem has been 
solved in the non-interacting case. 
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Our numerical results shows that the HTL screening mass, in the explored tem- 
perature range, turns out to be a few percent higher than the non-interacting 
value 27rT, (very) slowly approaching it from above as the temperature in- 
creases. 

We have also compared our results with the ones obtained in another analyt- 
ical approach based on a dimensional- reduced lagrangian [32,33], which also 
finds a positive correction of order g 2 T to the non- interacting result. Quite 
interestingly, once we expand perturbatively what comes out from our numer- 
ical findings for the asymptotic high-temperature behavior of the screening 
mass, we also get a correction to the free result of order g 2 T. The coefficient 
of this correction differs from the one found in [32,33] since, in the HTL ap- 
proximation, soft gluons corrections are neglected. 

Concerning the lattice data available so far, to our knowledge, all the screen- 
ing masses extracted from them approach the non-interacting value 2nT from 
below. This indeed appears in conflict with the studies performed in the contin- 
uum by us and, within a different framework, by the authors of Ref. [32,33]. 
Whether this is related to limitations of the present lattice results or by a 
breakdown of the approximations assumed in both the continuum studies 
(based on a separation of the scales T, gT and g 2 T, which is rigorously defined 
only in the regime g <C 1) is an open question. 

Indeed the positive correction found in the analytical approaches mainly arises 
from the thermal mass aquired by the quark modes at finite temperature, a 
fact which appears sound. On the other hand, the lattice results available so 
far do not extend up to very high temperatures. Hence it is possible that the 
apparent mismatch between the lattice and the continuum results is due to 
the fact that the weak coupling regime, which makes the approximations and 
the separation of momentum scales justified, is achieved only at temperatures 
not yet covered by the lattice studies. However a more careful investigation 
of the possible limitations of the lattice results available so far might be of 
interest. 
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A The running of the coupling 



Our numerical results refer to the case Nf = 2. The transition temperature, 
following [35], has been set to the value T c = 202 MeV. The ratio T c /A-^ = 
0.7721 for the Nf = 2 case is taken from Refs. [36,37,38]. The running of the 
gauge coupling is given by the two-loop perturbative beta-function, leading to 
the expression: 

g-\T) = 2b log - A- + £ log (2 log j (A.l) 

A WS °0 I. A MS J 

where 60 = j^s N-l — 2^^, b\ = p^pp ^102 — 38^^ The renormalization 

scale \x is usually taken to be of the order of the temperature, which, apart from 
Aqcd, represents the only physical scale entering into the problem. For what 
concerns its precise numerical value one should, in principle, let it vary within 
a reasonable range in order to get an estimate of the theoretical uncertainty of 
the calculation. Here, for the sake of simplicity, due to the huge numerical work 
required to produce our results, we adopt the choice \i = 1.1427rT suggested 
in [35]. 
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